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ABSTRACT 


A  satellite  has  been  observed  through  an  optical  telescope 
with  a  CCD  camera,  the  electrical  signal  stored,  and  later 
analyzed  with  the  aid  of  a  digital  computer.  The  unresolved 
satellite  image  is  contained  in  3  to  7  picture  elements  in  fifty 
contiguous  data  frames.  This  allows  a  better  understanding  of 
the  physical  characteristics  of  the  satellite  image  as  viewed 
through  the  atmosphere  by  an  electro-optical  detector. 

This  paper  describes  the  data  processing  programs  used  to 
analyze  the  data,  the  analysis  results,  and  implications  for 
signal  processing  design. 
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I .  INTRODUCTION 

A  CCD  camera  has  the  potential  of  improving  the  performance 
of  satellite  detection  systems  by  means  of  higher  quantum 
efficiency,  smaller  pixel  size,  and  more  pixels  on  a  focal  plane 
than  EBSICON  cameras.  An  efficient  satellite  detection  processor 
should  be  designed  to  exploit  the  signal  characteristics  of  the 
camera  it  uses  as  an  electro-optical  detector.  Selected  portions 
of  a  digital  satellite  detection  processor  (IMC  Signal  Processor) 
and  CCD  camera  are  being  developed  for  the  GEODSS  advanced  tech¬ 
nology  program.  The  output  of  a  CCD  camera,  when  observing  a 
satellite,  has  been  digitized,  and  then  analyzed  with  the  aid  of 
a  computer.  This  allowed  a  detailed  study  of  the  electrical  signal 
that  the  processor  would  have  to  use  prior  to  designing  the 
processor. 

Satellite  ATS-6  was  observed  with  the  31"  main  GEODSS  ETS 
telescope  at  the  White  Sands  Missile  Range.  The  CCD  camera  was 
at  the  prime  focus  and  was  operated  in  the  stare  mode  with  an  0.1 
second  integration  time.  The  ETS  photometer  measured  the 
satellite's  brightness  as  10.5  Mv.  The  camera  video  was 
digitized  to  12  bit  resolution,  sent  to  the  MODCOMP  computer 
at  the  ETS  and  stored  on  magnetic  tape.  This  tape  was  brought 
back  to  Lexington  and  placed  on  the  MODCOMP  computer  here.  The 
tape  was  read  by  the  MODCOMP  and  the  data  transmitted  to  the 
HP  9845  graphics  computer  and  stored  there.  The  HP  computer 
was  used  for  all  the  data  analysis. 


The  satellite  observation  is  contained  in  fifty  contiguous 
frames.  Each  frame  contains  an  array  of  32  x  32  picture  elements 
(pixels).  The  integration  time  for  each  frame  was  0.1  second. 
Each  pixel  has  the  equivalent  field-of-view  of  2"7  x  2’.’7. 
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II.  DATA  PROCESSING 

The  frames  of  data  were  read  from  the  digital  tape  by  the 

MODCOMP  computer  and  transmitted,  on  an  RS-232C  interface1,  to 

2 

the  HP  computer.  Using  this  interface  the  HP  computer  was 
operated  to  look  like  a  standard  terminal  to  the  MODCOMP.  The 
major  task  of  this  interface  was  to  write  an  HP  terminal 
emulator  program.  This  program  has  to  read  the  data  from  the 
interface,  store  the  data,  and  respond  to  the  MODCOMP  with 
standard  terminal  handshake  commands.  The  program,  "MOD-HP", 
listed  in  Appendix  A  did  this  at  9600  baud. 

To  operate  at  this  high  a  rate,  each  data  word  was  read 
as  a  character  string.  The  character  strings  were  stored  on 
HP  tape  cassettes  as  the  data  word  was  being  read.  The  1024 
character  strings  from  each  frame  were  separated  by  markers. 

Once  all  fifty  frames  of  the  character  strings  were  stored  on 
tape,  the  HP  computer  was  disconnected  from  the  MODCOMP,  and 
was  operated  independently  from  then  on. 

Two  tape  cassettes  were  required  to  store  all  the  character 
strings.  All  the  character  strings  were  then  read  off  the  tapes 
and  stored  on  one  floppy  disc  using  the  "COPY"  program  listed 
in  Appendix  A.  The  character  strings  cannot  be  used  in 
programs,  but  must  be  converted  to  the  numbers  they  represent. 
The  program  "READ",  Appendix  A,  read  the  character  strings  off 
the  disc,  converted  the  character  strings  to  numbers  and  stored 

IT!  Hewlett  Packard,  "Serial  I/O  Interface-98036A, "  1977. 

2.  Hewlett  Packard,  "I/O  ROM  Programming  System  45B,"  1979. 


these  on  tape.  The  "COPY"  program  was  again  used  to  transfer 
data  from  tape  to  disc. 

To  insure  that  the  data  stored  on  the  disc  was  accurate, 
a  print-out  of  the  data  stored  on  the  MODCOMP  tapes  was  compared 
to  the  data  stored  on  the  HP  disc.  The  data  was  identical  in 
frame- to- frame  and  entry-to-entry  comparisons.  The  first  data 
entry  in  each  frame  of  both  the  MODCOMP  and  HP  data  was  the 
same,  but  was  not  a  valid  data  value.  To  avoid  problems  in 
data  processing  the  value  of  the  first  data  entry,  in  each  frame, 
was  set  to  the  mean  of/  the  data  of  each  frame  and  is  then  called 
a  valid  data  frame. 

Now  the  data  is  ready  to  be  analyzed.  The  first  data 
analysis  program  is  "MAXAMP"  listed  in  Appendix  B.  "MAXAMP”  finds 
the  maximum  amplitude  in  each  frame  and  prints  the  frame  number, 
amplitude,  and  array  position.  A  sample  of  "MAXAMP"  output  is 
shown  in  Figure  1.  The  frame  number,  e.g.  T11V01:F  stands  for 
CCD- DIGITIZER  test  #11  (Til) ,  valid  data  frame  (V) ,  frame  number 
(01) ,  and  stored  on  floppy  disc  (F) .  The  pixel  array  position 
is  numbered  from  1  to  1024  in  each  data  frame. 

The  next  data  analysis  program  listed  in  Appendix  B  is 
"FREDIS"  which  plots  the  frequency  distribution  of  amplitudes 
for  each  frame.  Figure  2  contains  plots  from  two  of  the  fifty 
frames.  The  amplitude  scale  adjusts  for  each  frame  depending  on 
the  maximum  amplitude.  The  occurrence  scale  also  adjusts 


k 


FRAME 

T 1 1 V 10:  F 

MAXIMUM 

AMPLITUDE 

185 

POSITION 
2  7b 

T  1  1  V  1  1  :  F 

MAXIMUM 

AMPLITUDE 

182 

276 

T 11 V 1 2 : F 

MAXIMUM 

AMPLITUDE 

163 

27b 

t  1 1  v  1 3 :  F 

MAXIMUM 

AMPLITUDE 

159 

cl  i  i 

t  :  i  v  1 4 :  f 

MAXIMUM 

AM  PI.  I  TUDE 

224 

277 

T  1  1  V  1  5 :  F 

MAX  I  MUM 

AMPLITUDE 

176 

•*.  -?  -9 

C-  1  1 

T  1  1  V  1  6 :  F 

MAX  I  MUM 

AMPLITUDE 

288 

277 

i  n  v  1 7 :  f 

MAX  I  MUM 

AMPLITUDE 

173 

*7  **s 

k.  i  i 

T  1  1  V  1  3 :  F 

MAXIMUM 

AMPLITUDE 

147 

cL  f  i 

r  i  1  V  1  9 :  F 

MAXIMUM 

AMPLITUDE 

180 

d  i  r’ 

T 1 1 V20 :  F 

MAXIMUM 

AMPLITUDE 

172 

276 

T 1 l  V21 : F 

MAXIMUM 

AMPLITUDE 

143 

2  76 

T  1 1  V  2  2  :  F 

MAX  I  MUM 

AMPLITUDE 

144 

276 

r 1 1 v 23:  f 

MAXIMUM 

AMPLITUDE 

143 

2  {  6 

rt !  2  4 :  F 

MAXIMUM 

AMPLITUDE 

175 

2  76 

Ti  1  V  2  5 :  F 

MAXIMUM 

AMPLITUDE 

134 

276 

Fig.  1.  "MAXAMP"  sample  output# 


5 


depending  on  the  number  of  pixels,  excluding  the  zero  amplitude, 
for  which  there  are  from  350-400  pixels  in  each  frame. 

The  program  "PEAKS"  examines  a  frame  of  data  and  prints  out 
the  location  and  amplitude  of  any  pixel  whose  amplitude  is  above 
a  selected  level.  Figure  3  is  a  sample  of  "PEAKS"  output.  A 
program  listing  of  "PEAKS"  is  in  Appendix  B.  The  frame  number 
is  given  and  the  selected  level  (e.g.  19).  The  row,  column, 
array  position  (as  in  "MAXAMP"  output) ,  and  corresponding 
amplitude  are  listed. 

Program  "DATAST",  Appendix  B,  prints  out  the  amplitude 
values  of  any  frame  by  row  and  column.  An  asterisk  is  placed 
next  to  amplitudes  above  a  selected  amplitude.  The  output  is 
displayed  on  the  CRT  or  printed  on  the  line  printer.  A  sample  of 
output  from  frame  #11  is  shown  in  Figure  4.  The  top  line  of 
each  pair  is  array  elements  1  to  16  and  the  second  line  is 
array  elements  17  to  32.  In  this  frame  the  satellite  pixels 
have  amplitudes  of  182,  112,  32,  and  40.  An  isolated  pixel  of 
noise  with  amplitude  32  is  one  array  line  above  this. 

The  "FRAME"  program  (Appendix  B)  allows  a  13  x  11  pixel 
array  around  the  satellite  image  of  any  frame  to  be  displayed. 
After  the  array  of  pixels  is  displayed,  the  amplitudes  are 
plotted  first  in  column  and  then  in  row.  Figure  5  is  a  print¬ 
out  of  the  array  of  pixels  around  the  satellite  image  from 
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frame  #12.  Figure  6  is  a  plot  of  the  amplitude  in  the  sixth 
column  of  Figure  5.  The  dashed  line  is  a  level  at  an  amplitude 
of  19.  Figure  7  is  a  plot  of  amplitudes  in  the  fifth  row  of 
Figure  5.  The  dashed  line  is  also  at  an  amplitude  of  19. 

Appendix  C  contains  the  pixels  (5  x  6)  surrounding  the 
satellite  image  for  frames  1  through  12.  These  were  obtained 
using  a  modified  version  of  program  "FRAME" . 


III.  ANALYSIS  OF  DATA 

Analysis  of  the  output  of  the  "MAXAMP"  program  indicates 
that  the  highest  amplitude  of  the  satellite  image  moves  between 
two  adjacent  pixels.  Over  the  fifty  frames,  the  highest 
amplitude  switches  15  times  between  the  two  pixels,  spends  an 
equal  amount  of  time  in  each  and  there  is  no  periodicity  to  the 
switches.  The  range  of  highest  amplitudes  is  127  to  224  with 
a  mean  of  165.  The  mean  of  highest  amplitude  in  one  pixel  is 
162,  while  the  other  pixel's  highest  amplitude  mean  is  169. 

Examining  the  output  of  "FREDIS"  shows  that  in  all  frames 
most  of  the  pixels  have  amplitudes  from  0  to  17.  The  mean 
of  the  amplitudes  from  each  of  fifty  frames  ranges  from  4.4  to 
4.9.  This  is  strongly  due  to  the  large  number  of  pixels  (350  to 
400)  in  each  frame  with  amplitudes  of  zero.  The  distributions 
of  all  fifty  frames  are  similar.  All  frames  have  less  than  20 
pixels  above  an  amplitude  of  19. 

The  output  of  "PEAKS"  shows  that  in  every  frame  the  pixels 
with  amplitudes  greater  than  28  are  either  those  of  the  satellite 
image  or  of  isolated  pixels.  There  are  532  pixels  that  have 
amplitudes  greater  than  19  in  the  fifty  frames.  Most  of  these 
pixels  are  either  those  of  the  satellite  image  or  isolated 
pixels.  However,  12  of  the  pixels  in  10  of  the  fifty  frames 
have  an  adjacent  pixel  greater  than  19.  In  half  of  these  cases 
the  adjacent  pixels  are  vertically  related  and  in  half  horizontally. 
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The  pixels  with  amplitudes  greater  than  39  are  all  those  of 
the  satellite  image. 

The  output  of  "DATAST"  shows  that  the  noise  is  uncorrelated 
from  pixel  to  pixel  and  from  line  to  line.  Eight-four  percent 
of  the  pixels  with  an  amplitude  greater  than  8  have  no  adjacent 
pixel  with  amplitude  greater  than  8.  Ninety-three  percent  of 
the  pixels  with  amplitude  greater  than  12  have  no  adjacent  pixel 
with  amplitude  greater  than  12. 

The  output  of  "FRAME"  gives  information  primarily  about  the 
pixels  associated  with  the  satellite  image.  Table  1  is  a  summary 
of  the  distribution  of  satellite  image  pixels  by  frames.  The 
distribution  of  energy  from  the  satellite  image  is  very  dynamic. 
The  distribution  changes  significantly  in  periods  of  0.1  sec  to 
0.3  sec.  In  61.5%  of  the  frames  the  energy  distribution  doesn't 
fit  a  circular  Gaussian  function. 


IV.  SIGNAL  PROCESSING  CONSIDERATIONS 

After  reviewing  the  analysis  results,  it  is  apparent  that 
the  satellite  image  is  not  stable  over  time  periods  as  short 
as  0.2  sec.  The  instability  of  the  image  could  result  from  three 
primary  sources:  satellite  motion,  telescope  jitter,  and  atmospheric 
seeing.  During  the  time  of  data  collection,  the  ATS-6  satellite 
was  a  three-axis  stabilized  geosynchronous  satellite.  Its 
reflected  light  vs.  time  had  been  measured  by  several  observers 
at  the  ETS  and  was  found  to  be  steady  throughout  the  night,  as 
would  be  expected  from  this  type  of  satellite.  Therefore  satellite 
motion  didn't  account  for  the  image  instability  observed  in  the 
data.  The  fifty  frames  of  data  were  carefully  scrutinized  to 
detect  any  indication  of  telescope  jitter.  Telescope  jitter 
would  be  detected  in  the  data,  if  the  pixels  which  contain  the 
satellite  image  shifted  together  in  the  same  direction.  There  is 
no  observable  correlation  between  the  maximum  amplitude  pixel  shift 
direction  and  the  other  satellite  image  pixels'  shift  direction. 
Calculating  the  centroid  of  the  satellite  pixels  shows  that  the 
centroid  movement  is  in  random  directions  rather  than  in  one  or 
two  orthogonal  directions  that  would  characterize  telescope  jitter. 
Atmospheric  seeing  is  therefore  the  leading  candidate  as  the 
source  of  the  image  instability.  Although  few  measurements  of 
atmospheric  seeing  with  time  constants  as  short  as  0. 1-0.2  sec, 
and  with  two-dimensional  spatial  resolution  as  fine  as  2'.’7  have 
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3 

been  reported,  one  reference  does  report  the  characteristic 
image  fluctuation  seen  in  this  data  as  atmospheric  seeing.  Many 
other  reports  on  atmospheric  seeing  describe  image  fluctuation 
similar  to  that  observed  in  this  data.  Atmospheric  seeing  will 
affect  stellar  images  in  the  same  fashion.  A  signal  processing 
system  that  removes  stars  by  comparing  a  star  field  at  one  time 
with  the  same  star  field  at  a  time  0.2  sec  later  or  longer,  must 
incorporate  a  technique  to  reduce  false  alarms  due  to  the  star 
image  instability.  These  false  alarms  are  located  adjacent  to 
star  images,  and  a  guard  band  around  each  star  that  inhibits 
detections  would  remove  these  false  alarms. 

Although  the  resolution  of  the  telescope  is  1'.'0  or  smaller 
and  the  pixel  resolution  is  2'.'74,  the  satellite  image  is  in 
3  to  7  pixels.  Most  of  the  energy  is  in  2  to  3  pixels,  however 
in  74%  of  the  frames  the  energy  is  in  4  to  5  pixels.  In  more 
than  half  the  frames  the  distribution  of  satellite  energy 
doesn't  fit  into  a  circular  function.  There  is  little  correla¬ 
tion  from  frame- to- frame  of  the  spatial  distribution  of  satellite 
energy.  This  requires  that  predictions  of  satellite  SNR  and 
system  detection  performance  be  modified  to  incorporate  these 
physical  realities. 


3.  MIKESELL,  A.H.  et  al,  "The  Scintillation  of  Starlight",  JOSA, 
41,  689-695,  1951. 
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The  star  and  satellite  images  are  in  several  adjacent  pixels 
(at  a  2'.'74  pixel  size),  while  the  noise  is  only  one  pixel  long. 
This  is  an  excellent  discriminant  between  the  satellite  and  high 
amplitude  noise  that  competes  with  the  satellite  signal.  If 
a  threshold  level  of  28  was  used  with  the  pixels  in  the  fifty 
frames,  the  threshold  crossings  would  be  due  to  the  satellite 
(adjacent  threshold  crossings)  or  to  the  noise  (isolated  thres¬ 
hold  crossings).  The  isolated  threshold  crossings  could  then  be 
removed  and  only  the  satellite  threshold  crossings  would  remain. 
As  the  threshold  is  lowered,  noise  threshold  crossings  will  be 
more  frequent  and  false  alarms  due  to  noise  will  start  to  re¬ 
appear.  At  a  threshold  of  28,  100%  of  the  noise  false  alarms 
are  removed;  at  a  threshold  of  12,  93%;  and  at  a  threshold  of 
8,  84%.  The  remaining  false  alarms  due  to  noise  can  be  removed 
by  additional  processing  (streak,  slope,  velocity,  etc.). 
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V. 


CONCLUSIONS 


The  analysis  of  a  satellite  image  as  seen  by  a  CCD  camera 
mounted  on  the  31"  prime  focus  at  the  ETS  has  been  accomplished. 
The  complexity  of  the  image  as  a  function  of  spatial  distribution 
and  time  can  be  clearly  seen.  The  seeing  effects  were  worse  than 
what  had  previously  been  expected  at  the  ETS.  The  utility  of 
having  a  digitized,  stored  signal  from  a  sensor  before  designing 
a  signal  processor  that  uses  that  particular  sensor,  is  obvious. 

The  pixel  size  of  a  camera  is  important  to  the  detection 
system;  as  too  small  a  pixel  will  cause  the  satellite  energy  to 
be  spread  over  many  pixels  and  thus  reduce  the  SNR,  or  if  the 
pixel  is  too  big  the  night  sky  will  contribute  additional  noise 
and  also  lower  the  SNR.  The  image  size  depends  on  the  atmospheric 
seeing,  therefore,  the  selection  of  pixel  size  depends  on  the 
seeing  expected  at  a  site.  The  ideal  way  to  measure  seeing  at 
the  ETS  would  be  to  use  an  automated  version  of  what  was  done 
in  this  report.  The  microprocessor  based  PVDS  system  would  be 
ideal  for  this  purpose  and  could  provide  a  short  exposure 
measure  of  seeing.  Using  a  CCD  camera  mounted  on  the  cassegrain 
position  of  the  31"  ETS  telescope  would  provide  even  more  detail 
of  the  image  intensity  distribution.  Binary  stars  could  be 
viewed  to  eliminate  some  observing  uncertainties  and  data  could 
be  collected  on  a  regular  basis. 
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in 

I*  1  M  D»t  500  ■  [  362  1 

!  "MOD-HP" 

S'0 

'-RE  ATE  "TtlFOl", 1,10000 

30 

ASSIGN  #1  TO  "TllFOl" 

40 

RESET  10 

SO 

WRITE  10  10,5;  1 

60 

WAIT  WRITE  10, 4; 64 

ro 

WAIT  WRITE  10. 4; 203 

a  o 

WAIT  WRITE  10,4; 32+7 

90 

RE  AH  10  10,4;R4 

100 

WRITE  10  10, 5; O 

1  1  u 

WAIT  READ  10, 4; A 

1  READY 

1  2'0 

WRITE  10  lO. 7; 0 

i  TRIGGER 

1  30 

i 

1  4  0 

i  SET  UP  THE  9S0 

36A  TO  INTERRUPT  ON  INPUT 

DATA 

ISO 

! 

1  60 

WRITE  10  10,5; 128+4 

i  ro 

F  = I 0 F L A G •  lO  • 

1  30 

8=  1  OSTATUSt 10  • 

1  90 

PRINTER  IS  16 

0  O  O' 

PRINT  "10  FLAG  1  "  , F ,  "  I 

0  STATUS  1  “  ,  S 

2  1  0 

READ  10  10,5; R5 

090 

READ  1 0  1 0 , 6 ; R  6 

2  3  0 

STATUS  io; V 

24  0 

PRINT  “ R4  0  "  , R 4 ,  " R 5 

146  ",R5,"P6  248  ",R6, 

"  V 

4  02  ",  V 

2  S  0 

1  =0 

! INDEX  FOR 

WORDS 

COMING  OVER 

2  6  O 

DEEP 

270 

DEEP 

230 

TOPEN  10,2  GOSUE  Get  input 

2  3  0 

Loop: 

1  WAIT  FOR 

DAT  A 

FROM  MO DC  CMP 

3  O'  0 

ON  LED  1  GOTO  310 

3  1  0 

L  *  =  L  B  D  * 

320 

E  X « P 0 S  <  K  $ ,  c HR*  1  255  >  &CHR*  •;  2 

1  >  > 

3  3  0 

IF  E  •.  0  THEN  440 

3  4  0 

GOTO  Loop 

350 

Ge  t i nput  ; 

!  GET  INPUT 

FROM 

INTERFACE 

360 

1  =  1  +  1 

370 

IF  I  500  THEN  Loop 

380 

D  $ ( I  > =TEUF * 

3  9  0 

I  F  PC'S  <.  D *  (.  I  :>  ,  CHR*  t.  5  >  >  THEN 

OUTPUT  10; CHR* < 6 > 

4  0  0 

PRINTER  IS  16 

410 

PRINT  I 

420 

WAIT  900 

4  30 

RETURN 

440 

FOR  1=1  TO  50O 

450 

PRINT  *1 ; I , D*< I > 

460 

NEXT  I 

470 

END 

"MOD-HP" 
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APPENDIX  A 


FOR  k-J  =  1  0  TO  50  !  "i?EAD" 

D  I  M  Df  (.  2O0  >  C  320]  ,  PC 321  > ,  0<  200 >  .  F  <  1  024  > 

X*=VAL*‘ W> 

M 1 “ T 1 1F"8.X*« " : F8 " 

M2*  =  "T1 1D"&X$«,"  :  T15" 

FISSION  #1  TO  Wl* 

ON  END  #1  GOTO  680 
St  =0 
D  =  0 

FOR  1-1  TO  20O  !  READ  STRING  FROM  TAPE  cl' 

REFID  #i;  I,D*<  I  > 

D*  < ! >«D# C I > C 1 , 31 93&CHR*< 1 3 > 

PRINTER  IS  16 

FOR  ,T=  1  TO  200  !  BREAK  STRING  INTO  VARIABLE 

IF  j: 1  THEN  290 
F  <  1  =POS  <  B*<  I  >  ,  CHR*  <1  3  >  ) 

L f  =  P 0 S ( D* < I >  C 1 , P < 1)3, CHR$ < 1 0> ) 

L  =  LEN>  D*C I  * C 1 , PC  1)  3 > 

IF  Lf<>0  THEN  240 

IF  L<2  THEN  280 

OC  1  >=VALCD*C  I  )U  ,  PCI  )  3  > 

D  =  1 

GOTO  290 

IF  L<  3  THEN  280 

0C1)=VALCD*CI>C2,PCD3> 

D=1 

GOTO  290 
0  c 1 )=-l 
Jl= J+l 

P  <  .T 1  )  =POS  <  D*  C  I )  [  P  <  J  >  + 1  3  ,  CHR*<  1 3)  ) 

IF  P C J 1 ) =  0  THEN  480 
PCJ1>=PCJ1 >+PCJ> 

L=LENC  D$  C  I  >CPC J>  +  1 , PC  J1)  3  > 

IF  L < 4  THEN  460 

Lf  =  POS C  D* CI>, CHR$  < 10 ) > 

IF  Lf < : 0  THEN  390 
G=  1 

GOTO  430 
G  =  2 

Dt<  I  )  t  P  C  J  >+G,  PC  J 1  >  3  =  TRI  M  $  <  D  tU  )  [  F'  <  J  >  +G ,  P  C  J  1)3) 


APPENDIX  A 


4  id 
4  2  0 
430 
440 
450 
460 

4  70 
480 
490 
500 
510 
520 
530 
540 
550 
560 

5  70 
580 
590 
600 

6  1  0 
620 
6  3  0 
640 

6  5  0 
660 
670 
680 
690 
700 
710 
720 

7  30 
740 
750 
760 
778 
780 
790 
80O 


L=LEN<  D#< I>[P<  J>  +  G,P< Jl>  ]  :» 

IF  L  6  THEN  460 

0  <  J 1 >  * VHL  <  D*  < I ) C  P  <  J  >  +C , P  <  J 1 >  3  ) 

Ii  =  H  -f  1 

GOTO  470 

0  <  J 1  '=-1 

GOTO  500 

S  =  D 

GOTO  510 
NEXT  J 


Ss  =  0 

C**St + 1 -Si 
C  f  =  3 1 +8 


n=e 

IF  < 0 <  1  >  =  - 1  >  AND  <S=0>  THEN  670 
F U R  K  —  C i  TO  C f 
L  =k-St  -t-Ss 

IF  0<L>=-1  THEN  608 

GOTO  640 

L  =  L  +  1 

Cf  =  C  f -1 

S  s  =  8  s  ♦  1 

GOTO  588 

F  K  :  •  =  0  <  L  > 

NEXT  K 
St =  St  +  S 
NEXT  I 

FOR  f'1=  1  TO  32 

FOR  N=t1  TO  102  3+M  STEP  32 

I  URGE  #,4D 

PRINT  USING  700;FCN> 

NEXT  N 

PRINT  L I N < 1 > 

NEXT  M 


ON  END  #2  GOTO  790 
CREFITE  W2*,l,,9000 
ASSIGN  #2  TO  142$ 
PRINT  #2JF(*> 

NEXT  W 
END 


"READ"  (Cont.) 


NEXT  J 


PACK  VARIABLES  FOR  OUTPUT 


NEXT  K 
NEXT  1 

OUTPUT  DATA  <M> 


NEXT  M 


10 

FOR  1=10  TO  50 

!  "COPY" 

20 

I*=VAL#<  I  ) 

!  FRAME  NUMBER 

30 

T$="T1 ID" 

40 

Pl*=" : T15" 

!  TAPE  DRIVE 

50 

P2#=" : F" 

!  DISC 

60 

F1$=T*^I$&P1* 

70 

F2#  =  T#U*&P2# 

80 

PRINT  F1#,F2# 

90 

COPY  FI#  TO  F2# 

100 

NEXT  I 

110 

END 

"COPY" 
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10 

PRINT  "FRAME","  ","  “."POSITION" 

!  "MRXAMP"  ! 

20 

FOR  1=10  TO  50 

30 

DIM  FO024) 

!' 

40 

I #  =  VAL#  < I ) 

50 

u_ 

* 

M 

> 

W 

b- 

II 

** 

60 

ASSIGN  #1  TO  V# 

70 

ON  END  #1  GOTO  100 

80 

READ  *1;F<*> 

90 

Ma>:=0 

i 

100 

FOR  M= 1  TO  32 

!  OUTPUT  DATA  <M>  i 

110 

FOR  N=M  TO  1023+M  STEP  32 

1 

120 

IF  F < N ) >Max  THEN  140 

130 

GOTO  160 

140 

Max=F  <N) 

150 

Nax=N 

( 

1  60 

NEXT  N 

1 

170 

NEXT  M 

NEXT  M 

1  80 

PRINT  VS, "MAXIMUM  AMPL I  TUBE " , Max , Nax 

i 

190 

PRINT  LINtl) 

I 

200 

NEXT  I 

i 

210 

END 

1 

"MAPCAMP" 


10  PRINTER  IS  16  !  "FRED IS" 

20  PRINT  PAGE;"  This  program  allows  the  user  to  enter  a  frame  of  data  and  ha 
"e  a" 

3©  PRINT  "  frequency  distribution  computed  and  drawn." 

43  PRINT  LIN<1>;"  After  the  distribution,  a  beep  will  sound.  The  chart  will 


5£i  PRINT  "remain  on  the  screen  and  the  program  will  halt  until  you  press  CONT 

M 

60  PRINT  "This  will  give  you  time  to  digest  the  information  on  the  screen." 

70  PRINT  LINO);"  If  you  want  the  output  of  this  program  dumped  from  the  CRT 

II 


S3  PRINT  "to  the  internal  thermal  printer,  then  answer  Y  to  the  question  below 


90  PRINT  "If  you  don't  want  hardcopy  answer  N.  Then  press  CONT." 

100  INPUT  "Do  you  want  hardcopy  <YxN>?" , At 

110  IF  UPC#<A#)="Y"  THEN  Yhard 

120  IF  UPC# <  At )  =  " N "  THEN  Nhard 

13©  BEEP 

140  GOTO  100 

150  Yhard:  Hardcopy=l 

160  GOTO  180 

170  Nhard:  Hardcopy=0 

130  N=224 

190  CALL  Dr i ver<N, Hardcopy > 

200  END 

210  SUB  Dr i ver < N, Hardc opy ) 

220  RANDOMIZE 

2  30  DIM  A<N),L#<N),Units#t30],V<1024) 

"FREDIS" 
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24©  INPUT  "Which  frame  01  to  50?“,  Fr* 

250  V*  =  " T 1 1 V " &Fr  " : F" 

260  ASSIGN  #1  TO  V* 

270  READ  #1  ;v<*> 

2 SO  Ma**V< 1 > 

290  FOR  1=1  TO  1024 
SO0  Max=MAX *  Max , V ( J  > > 

510  NEXT  J 

520  MAT  A=ZER 

550  H i gh  =  0 

540  Lh i gh=4096 

350  FOR  K = 1  TO  1024 

560  IF  i  V  f  K  >  >H  l  gh  >  AND  <  V  *  k  >  L h  i  g h  )  THEN  380 
3?0  GOTO  390 
3 SO  H  i  g h = V  (.  K  > 

590  NEXT  K 

40O  FOR  J 1 = 1  TO  1024 

418  IF  V1'  J  1  >  =  H  i  gh  THEN  43© 

420  GOTO  44© 

430  A  <  H  i  gh  )  =A'.Hi  gh  >  +  1 

440  NEXT " J 1 

450  IF  High*®  THEN  490 

460  Lhigh=High 

470  High=0 

480  GOTO  350 

496  PRINTER  IS  0 

500  PRINT 

510  PRINT  V* 

520  PRINT  "MAX  AMP" 

530  PRINT  Max 
540  PRINTER  IS  16 

550  FOR  1=20  TO  Max  STEP  20  !  ENTER  THE  TITLES 

560  L  $ < I > =  VAL* ( I > 

570  NEXT  I 

5 8 O  IJ n its**" N u nt b e r  o f  Pixels" 

59©  T i 1 1 *#* " Amp 1 i t  udes " 

6  O  O  N  =  M  ax 

610  CALL  B ar <  A  <  *  )  ,  L  * < *),N, Units*, Title* > 

620  BEEP 
6  3  O  P  A  IJ  S  E 

640  IF  Hardcopy* 1  THEN  DUMP  GRAPHICS 
650  EXIT  GRAPHICS 
6  6  O  S  U  E  E  X I T 

670  SUB  B  ar  <1 A  ( *  )  ,  L*  <  *  >  ,  N ,  Un  i  t  s* ,  T  i  ♦  1  e*  ) 

680  PLOTTER  IS  "GRAPHICS" 

690  GCLEAR 
70©  GRAPHICS 
710  PEN  1 
720  Sum*© 


"FREDIS"  (Cont.) 
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7  50 

Max=Rt 1 > 

740 

FOR  1=1  TO  N 

750 

Sum=Sum+A<  1  > 

760 

Ns*:  =  MRX < Mix ,  fl<  I  >  ) 

770 

NEXT  I 

780 

Ly=LGT  <•  Max  ) 

780 

SETGU 

300 

LOCATE  15,120,40,95 

!  Set  screen  boundaries  and  scale 

810 

CLIP  0,123,40,100 

850 

MOVE  61.5,97 

830 

LORG  5 

840 

LABEL  USING  850; Ti t 1 e* 

850 

IMAGE  # , K 

860 

SCALE  0,N+.5,0, M ax + Max ' 1 0 

370 

Xt i c  =0 

!  No  t  i  c 

mark s  on  hoc i zont al  ax  i  s 

880 

Test,  y  t  i  c  =FRACT  <  Ly  )  +  < Ly<  0  > 

890 

Yt  i  c  =  10' ( I NT (Ly  > - 1 ) *< 1  + 1 . 5* 

<<Test.ytic>. 39794)  AND 

<  T est  yt  i  c  <.  .  69897  )  )  +4*  <  <  T e 

St  yt 

ic>=. 69897)  AND  < Test yt i c < =. 

87506) >+6.5*<Testytic> 

.  87506)) 

900 

!  Compute 

t  he  vert i  c  al  tic  mark  s 

910 

CALL  Laxes < Xt i c , Yt i c , 0, 0 ,  1 , 

1,2,0,  N-t.  5,0,  Max+Maxx 1 0 > 

920 

LORG  4 

930 

FIXED  1 

940 

FOR  1=0  TO  N 

950 

MOVE  I - . 5, 0 

960 

I PLOT  0 , A < I > , - 1 

!  Draw  the 

bar  for  the  ith  data  item 

970 

I PLOT  1,0 

980 

I PLOT  0 ,  -  A  d  ) 

990 

MOVE  I, A(I) 

1000 

NEXT  I 

1018 

DEG 

1020 

IMAGE  K 

10  30 

LDIR  90 

!  Rotate 

the  label  90  degrees 

1040 

SETGU 

1050 

LORG  5 

1060 

MOVE  5,70 

1070 

LABEL  USING  1020; Units* 

!  Write 

the  display  units  being 

1080 

!  used  along  the  vertical  axis 

1090 

CLIP  0,123,0,100 

1100 

SETUU 

1110 

FUR  1=1  TO  N 

"FREDIS"  (Cont.) 
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f 


1120  MOVE  I, -Max.  20 
1130  LORG  8 
1 140  PEN  1 

1150  LDIR  90  !  Rotate  the  labels  98  degree 

1180  LABEL  USING  1020; L*<  I) 

1170  NEXT  I 
1180  BEEP 
1190  SUBEXIT 

1200  SUB  Laxes<Xt-  i  c  ,  Yt  i  c  ,  Xorg,  Yorg,  Xmaj  ,  Ymaj  ,  M  i  rit  i  c  s  i  ze ,  Xm  i  n,  Xmax ,  Ym  i  n ,  Ymax 
1210  IF  (Xmi n>*Xmax)  OR  < Ymi  n>=Ymax>  THEN  SUBEXIT 
1 220  Xf udge  = . 02*  <  Xmax-Xm i  n ) 

1 230  Yf udge= . 02* < Ymax-Y min) 

1 240  AXES  Xt  i  c  ,  Yt  i  c  ,  Xorg ,  Yorg ,  Xmaj  ,  Ymaj  ,  M  i  nt  i  c  s  i  ze 
1250  DEG 

1260  IF  NOT  Xt i c  THEN  Labely 

1270  LabelxlIF  SGN<Xtic)*-l  THEN  Parx 

1280  LDIR  90 

1290  LORG  8 

1300  GOTO  1330 

1310  Parx:  LDIR  0 

1320  LORG  6 

1330  FOR  I=Xorg  TO  Xmax  STEP  ABS(Xtie) 

1340  MOVE  I , Yor g-Yf udge 
1350  LABEL  USING  1360; I 
1360  IMAGE  #,K 
1370  NEXT  I 

1380  Lately:  IF  NOT  Yt i c  THEN  SUBEXIT 

1390  IF  SGN<Ytic)*-l  THEN  Pary 

1400  LDIR  0 

1410  LORG  8 

1420  GOTO  1450 

1430  Pary!  LDIR  -90 

1440  LORG  6 

1450  FOR  I = Yorg  TO  Ymax  STEP  ABS< Yt i c  > 

1460  MOVE  Xorg-Xfudge, I 
1470  LABEL  USING  1360;  I 
1480  NEXT  I 
1490  SUBEXIT 


"FREDIS"  (Cont.) 
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10 

DIM  F Cl 024) 

!  "PEAKS" 

20 

F$= " T 1 1 V22 : F  " 

30 

L=  1 9 

40 

PRINTER  IS  0 

50 

PRINT  FJ ,  "PEAKS  ABOVE " , L 

60 

PRINT 

70 

PRINT  "ROW", "COLUMN” , "POSITION". 

"AMPLITUDE" 

80 

PRINTER  IS  16 

80 

ASSIGN  #1  TO  F$ 

100 

ON  END  #1  GOTO  120 

1  10 

READ  #i;F<*> 

120 

FOR  M= 1  TO  32 

!  OUTPUT  DATA 

1  30 

C  =  0 

140 

FOR  N=M  TO  1023+M  STEP  32 

150 

C=C  + 1 

160 

IF  F  <  N  )  >  L  THEN  200 

170 

IMAGE  #, 4D, "  " 

180 

PRINT  USING  170;FCN) 

190 

GOTO  250 

200 

IMAGE  # ,  4  D ,  “  *  " 

210 

PRINT  USING  200;FCN) 

220 

PRINTER  IS  0 

230 

PRINT  M,C,N,F<N> 

240 

PRINTER  IS  16 

250 

NEXT  N 

260 

PRINT  LINC1) 

270 

NEXT  M 

!  NEXT  M 

280 

PRINTER  IS  0 

290 

PRINT  LINO) 

300 

END 

"PEAKS" 
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10  DIM  F < 1 024 )  !  " DAT AST " 

20  INPUT  11  AMPL  I  TUBES  HIGHER  THEN"  ,  L 
30  FOR  J=10  TO  50 

40  PRINT  LINO) 

50  PRINT  "FRAME  NUMBER ",J 

60  PRINT  LINO) 

70  S*=VAL*OJ> 

80  ASSIGN  #1  TO  "TllVlSn  " !  F" 

90  ON  END  #1  GOTO  110 

100  READ  #i;FO) 

110  FOR  M=1  TO  32  !  OUTPUT  DATA  <M> 

120  FOR  N=M  TO  1023+M  STEP  32 

130  IF  F  <  N >  > L  THEN  170 

140  IMAGE  #,4B,"  " 

156  PRINT  USING  I 40; F<N> 

160  GOTO  190 

170  IMAGE  #, 4D, 

180  PRINT  USING  170JFCN) 

198  NEXT  N 

200  PRINT  LINO  ) 

210  NEXT  M  !  NEXT  M 

220  NEXT  J 

230  END 


"DATAST" 
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10  INPUT  “LEVEL " , L  !  “FRAME" 

20  FOR  1=10  TO  50 

30  DIM  F  <  1 024  > 

40  I*=VAL*U> 

50  V  $  = " T 1 1 V " & I ♦ & " : F " 

60  ASSIGN  II  TO  Vt  !  DATA  FROM  DISC 

70  ON  END  #1  GOTO  90 

80  READ  # 1 ; F  <  *  > 

90  GRAPHICS 

100  LINE  TYPE  1 

HO  SCALE  0,12,0,15 

128  FOR  M= 1 5  TO  25  !  LABEL  AMPLITUDE 

130  C=0 

140  FOR  N=M  TO  1023+M  STEP  32 
150  C=C+1 

160  IF  FCNKL  THEN  200 
170  MOVE  M- 1 5 , C 

180  LABEL  F(N) 

190  IF  C=  1  3  THEN  210 

200  NEXT  N 

210  NEXT  M 

220  WAIT  3000 

230  GCLEAR 

240  SCALE  0, 13,0,230 

250  FOR  M= 19  TO  22  !  PLOT  COLUMNS 

260  WAIT  200 

270  C=0 

280  GCLEAR 

290  MOVE  0 , L  !  PLOT  LEVEL 

300  LINE  TYPE  4 

310  DRAW  13, L 


APPENDIX  B 


320  LINE  TYPE  1 

330  MOVE  1,150 

340  LABEL  “COLUMN" 

350  MOVE  0.0 

360  FOR  N*M  TO  1023+M  STEP  32 

3?0  C  =  C+  1 

380  PLOT  C-l,F<N> 

390  PLOT  C,F< N) 

400  NEXT  N 

410  WAIT  500 

420  MOVE  0,0 

430  NEXT  M 

440  GCLEAR  ,  PL0T  ROWS 

450  MOVE  0,0 

4S0  FOR  N 1=239  TO  335  STEP  32  ,  rLOt  LEVEL 

470  MOVE  0 ,  L 

480  LINE  TYPE  6 

4<10  DRAW  1  3  ,  L 

500  LINE  TYPE  1 

510  MOVE  1,150 

520  LABEL  “ROW" 

530  MOVE  0,0 
540  FOR  M* 1  TO  13 
550  N=N1+M 
560  PLOT  M- 1 , F < N > 

570  PLOT  M,F<N> 

580  NEXT  M 

590  MOVE  0,0 

600  WAIT  500 

610  GCLEAR 
620  NEXT  N 1 

630  NEXT  I 

64W  END 


"FRAME"  (Cont.) 
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